AD-A219  353 


NUSC  Technical  Document  8589 
21  November  1989 


Some  Operational  Formulations 
For  Computational  Solutions 
Of  Underwater  Acoustic 
Scattering  Problems 


Environmental  and  Tactical  Support  Systems  Department 


DTIC 


Naval  Underwater  Systems  Center 

Newport,  Rhode  Island/New  London,  Connecticut 


Approved  for  public  release;  distribution  is  unlimited. 

no  on  is  178 


Preface 


This  document  was  prepared  for  the  Office  of  the  Associate  Technical  Director  for 
Technology  by  Dr.  David  Middleton. 


Reviewed  and  Approved:  21  November  1989 


W.  A.  Von  Winkle 

Head,  Environmental  and  Tactical  Support  Systems  Department 


REPORT  DOCUMENTATION  PAGE 


14.  REPORT  SECURITY  CLASSIFICATION 


1b.  RESTRICTIVE  MARKINGS 


24.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OECLA&SIFICATION/  DOWNGRADING  SCHEDULE 


2.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  public  release; 

distribution  unlimited 


4.  PERFORMING  ORGANI2ATION  REPORT  NUMBER(S) 
TD  8589 


S.  MONITORING  ORGANIZATION  REPORT  NUMBERCS) 


6*.  NAME  OF  PERFORMING  ORGANIZATION 
Naval  Underwater  Systems  Cent4 


6b  OFFICE  SYMBOL 
(If  * pplktbl*) 

r  Code  31 


74.  NAME  OF  MONITORING  ORGANIZATION 


Be  ADDRESS  (C/fy,  Stitt,  tnd  ZIP  Cod*) . 

New  London  Laboratory 
New  London,  Conn.  06320 


7b.  ADDRESS  (Gfy.  State.  md  ZIP  Cod*) 


Be.  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL 
Of  ipplkibl* ) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


8c.  ADDRESS  (City,  Sttf.ind  ZIP  Cod*) 


10  SOURCE  OF  FUNOING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

TASK 

WORK  UNIT 

NO. 

NO. 

ACCESSION  NO 

1 1  TITLE  (Include  Security  Oiaifkitton ) 

Some  Operational  Formulations  for  Computational  Solutions  of  Underwater  Acoustic 
_  Scattering  Problems _ 


12.  PERSONAL  AUTHOR(S) 

David  Middleton 


13a.  TYPE  OF  REPORT 

Tech.  Document 


13b.  TIME  COVERED 
FROM  9^88  TO  8/89 


16.  SUPPLEMENTARY  NOTATION 

See  Glossary  of  Principal  Symbols 


114.  DATE  OF  REPORT  (Year,  Month.  Oiy) 
1989  November  21 


IS  PAGE^l 


:ount 


17 

COSATI  CODES  S 

FIELD 

GROUP 

SUB-GROUP 

14  ABSTRACT 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  nvmbtr) 

^.Operational  Formulations  for  Propagation  and  Scattering; 
Underwater  Acoustic  Scattering  Problems;  Inhomogeneous 
determinigt-ic  and  random  media;  Langevin,  Dyson,  Bethe- 


Scattering  from  ocean  wave  surfaces,  volumes,  and  bottoms  remains  a  perennial  problem 
not  only  in  underwater  acoustics  but  also  in  the  use  of  acoustical  methods  in  remote  sensing  and 
modeling  of  the  medium  itself.  Although  analytic  solutions  are  sometimes  available,  cf.  [1],  [6],  in 
many  important  cases,  particularly  those  involving  multiple  scatter  and  non-isogradiant  (Vc  *  0) 
phenomena  including  the  role  of  both  stationary  and  moving  arrays,  and  general  geometries, 
effective  modeling  and  quantitative,  predictive  results  ultimately  require  computational  methods. 
This  is  true  for  both  deterministic  and  random  media. 

It  is  the  aim  of  this  introductory  study  to  present  a  combined  analytical -operational 
framework  from  which  computational  solutions  to  underwater  acoustic  scattering  problems, 
including  relevant  ocean  boundary  interactions,  can  be  initiated  and  carried  out.  It  is  emphasized 
that  here  explicit  computational  procedures  are  not  developed.  Rather,  the  specific  aim  is  to  present 


(cont'd) 


(cont'd) 


20  distribution /availability  of  abstract 

□  UNCLASSIFIE D/UNLIMITED  □  SAME  AS  RPT 


21. 


ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


00  FORM  1473. 84  mar 


81  APR  edition  mey  be  used  until  exhausted. 
All  other  editions  ere  obsolete. 


SECURITY  CLASSIFICATION  op  this  page 
Unclassified 


lecuwirv  CLASiiriCATie*  or  this  p«ck 
^  Continuation  of  Block  18: 


continuation  ox  ±uocx  ,,  ■  ■ 

Salpeter  Equations;  perturbation  series  solutions;  Born,  Bourret  /approx,  ;  Formulations 
for  computational  solutions;  feedback  operational  solutions.- 

Continuation  of  Block  19: 

the  "macro- algorithms"  from  and  to  which  appropriate  software  may  be  designed  and  applied. 

This  is  done  by  providing  formal  operational  structures  (with  some  special  cases  in  which  analytic 
solutions  are  obtainable).  A  collateral  aim  of  this  work  is  to  entice  interest  in  pursuing 
computational  formulations,  based  on  the  various  formal  operational  "solutions"  derived  here  for 
both  deterministic  and  general  random  acoustical  (i.e.,  scalar)  fields,  a(R,t). 

For  inhomogeneous  deterministic  media  (Vc  *  0,  c  «  c(r,t)  =  c©  +  £(r,t),  with  t 
nonrandom,  for  example),  operational  solutions  may  be  derived  directly:  all  representations  of  the 
medium  (and  boundaries)  are  equivalent.  The  resulting  solutions  provide  a  possible  alternative 
approach  to  the  current  "standard"  methods  involving  partial  or  parabolic  equations  (PE) 
developments  [7],  [8],  at  the  same  time  including  general  arrays,  geometries,  and  signals. 

On  the  other  hand,  most  of  the  scattering  mechanisms  in  the  ocean  are  essentially  random. 
Accordingly,  in  these  cases,  scattering  "solutions"  are  based  on  the  appropriate  Langevin  equation 
governing  propagation,  namely,  on  the  ensemble  of  representations  (a(R,t))  for  which  various 
appropriate  probability  measures  (often  nongaussian)  are  developed,  based  on  both  the  physics, 
geometry,  and  the  random  mechanisms  involved.  Thus,  practical  solutions  here  are  the  various 
moments  coo,  «xi0C2>,  etc.,  and  pdfs  wj(a),  w^a^ctj),  etc.,  which  are  needed,  for  example, 
in  the  calculation  of  scattering  cross  sections  and  subsequent  signal  processing  for  target  detection 
and  remote  sensing  [2]— £4], 

The  desired  framework  is  obtained  by  recognizing  that  both  the  deterministic  and  random 
scattered  fields  in  these  linear  media  are  representable  as  four-dimensional  feedback  systems,  with 
the  inhomogeneous  portions  (including  the  ocean  interfaces)  constituting  the  feedback  operations. 

New  features  of  this  work  include  a  combination  of  diagram  methods,  pertuibational  operational 
series,  and  equivalent  (i.e.,  deterministic)  media  operators,  which  then  provide  the  needed 
structure  for  both  the  direct  and  moment  solutions  a,  or  <aj>,  etc.,  independently,  with  the  latter 
based  on  the  resultant  Dyson  and  Bethe-Salpetcr  forms,  for  <a>,  <ajCi2>,  etc.[2],  [3].  Various 
analytical,  iterative,  and  approximative  examples  and  procedures  are  briefly  presented,  to  illustrate 
the  desired  formal  solutions  and  to  provide  elements  of  a  framework  for  computational 
implementation. 

The  interfacing  between  these  formally  analytic  solutions,  their  software  realizations,  and 
the  acquisition  of  (nearly)  on-line  numerical  results  remain  a  major  challenge. 

Aoeasslea  far _ 

~WIS  *  ORAAX 
DTIC  TAB  TJ 

UnannounoeA  □ 

Justlfleatlon _ _ 


9f - 

(  Plat rl but ion/ _ 

Availability  Codes 
*”  ~  Avail  and/or- 
91st  Speolal 


iccuntTv  CLAftsmcATie*  or  this  p»gc 


TD  8589 


Some  Operational  Formulations  for  Computational  Solutions  of 
Underwater  Acoustic  Scattering  Problems* 

by 

David  Middleton** 

Abstract 

'''  Scattering  from  ocean  wave  surfaces,  volumes,  and  bottoms  remains  a  perennial  problem 
not  only  in  underwater  acoustics  but  also  in  the  use  of  acoustical  methods  ir.  remote  sensing  and 
modeling  of  the  medium  itself.  Although  analytic  solutions  are  sometimes  available,  cf.  [1],  [6],  in 
many  important  cases,  particularly  those  involving  multiple  scatter  and  non-isogradiant  (Vc  *  0) 
phenomena  including  the  role  of  both  stationary  and  moving  arrays,  and  general  geometries, 
effective  modeling  and  quantitative,  predictive  results  ultimately  require  computational  methods. 
This  is  true  for  both  deterministic  and  random  media. 

It  is  the  aim  of  this  introductory  study  to  present  a  combined  analytical-operational 
framework  from  which  computational  solutions  to  underwater  acoustic  scattering  problems,  . 
including  relevant  ocean  boundary  interactions,  can  be  initiated  and  carried  out.  It  is  emphasized 
that  here  explicit  computational  procedures  are  not  developed.  Rather,  the  specific  aim  is  to  present 
the  "macro-algorithms"  from  and  to  which  appropriate  software  may  be  designed  and  applied.  <  t 
This  is  done  by  providing  formal  operational  structures  (with  some  special  cases  in  which  analytic 
solutions  are  obtainable).  A  collateral  aim  of  this  work  is  to  entice  interest  in  pursuing 
computational  formulations,  based  on  the  various  formal  operational  "solutions"  derived  here  for 
both  deterministic  and  general  random  acoustical  (i.e.,  scalar)  fields,  cc(R,t). 

For  inhomogeneous  deterministic  media  (Vc  *  0;  c  =  c(r,t)  =  Co  +  £(r,t),  with  c 
nonrandom,  for  example),  operational  solutions  may  be  derived  directly:  all  representations  of  the 
medium  (and  boundaries)  are  equivalent.  The  resulting  solutions  provide  a  possible  alternative 
approach  to  the  current  "standard"  methods  involving  partial  or  parabolic  equations  (PE) 
developments  [7],  [8],  at  the  same  time  including  general  arrays,  geometries,  and  signals. 

On  the  other  hand,  most  of  the  scattering  mechanisms  in  the  ocean  are  essentially  random. 
Accordingly,  in  these  cases,  scattering  "solutions"  are  based  on  the  appropriate  Langevin  equation 
governing  propagation,  namely,  on  the  ensemble  of  representations  (a(R,t)}  for  which  various 
appropriate  probability  measures  (often  nongaussian)  are  developed,  based  on  both  the  physics, 

*  Based  in  part  on  earlier  work  supported  by  the  Office  of  Naval  Research  (1970-1982)  and  on  more  recent  studies  for 
Codes  10  and  31,  NUSC,  1988,  1989.  See  also  [25]  herein. 

**  Contractor,  Physics  and  Applied  Mathematics,  127  E.  91  Street,  New  York,  NY  10128,  USA. 
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geometry,  and  the  random  mechanisms  involved.  Thus,  practical  solutions  here  are  the  various 
moments  «x>,  <ai<X2>,  etc.,  and  pdfs  wj(a),  wi(ai,02),  etc.,  which  are  needed,  for  example, 
in  the  calculation  of  scattering  cross  sections  and  subsequent  signal  processing  for  target  detection 
and  remote  sensing  [2]— {4]. 

The  desired  framework  is  obtained  by  recognizing  that  both  the  deterministic  and  random 
scattered  fields  in  these  linear  media  are  representable  as  four-dimensional  feedback  systems,  with 
the  inhomogeneous  portions  (including  the  ocean  interfaces)  constituting  the  feedback  operations. 
New  features  of  this  work  include  a  combination  of  diagram  methods,  perturbation al  operational 
series,  and  equivalent  (i.e.,  deterministic)  media  operators,  which  then  provide  the  needed 
structure  for  both  the  direct  and  moment  solutions  a,  or  <aj>,  etc.,  independently,  with  the  latter 
based  on  the  resultant  Dyson  and  Bethe-Salpeter  forms,  for  <a>,  «Xi<X2>,  etc.[2],  [3].  Various 
analytical,  iterative,  and  approximative  examples  and  procedures  are  briefly  presented,  to  illustrate 
the  desired  formal  solutions  and  to  provide  elements  of  a  framework  for  computational 
implementation. 

The  interfacing  between  these  formally  analytic  solutions,  their  software  realizations,  and 
the  acquisition  of  (nearly)  on-line  numerical  results  remain  a  major  challenge. 


1.  Introduction 

Scattering  from  ocean  wave  surfaces,  volumes,  and  bottoms  is  an  enduring,  central 
problem  not  only  in  underwater  acoustics  but  also  in  the  use  of  acoustical  methods  in  remote 
sensing  and  modeling  of  the  medium  itself.  Fortunately,  in  a  number  of  important  applications 
(where  first-order  Bom  approximations  are  adequate),  analytic  solutions  arc  available  (see,  for 
example,  recent  work  of  the  author  [1],  [4]). 

In  many  cases,  however,  strictly  analytic  methods  fail,  primarily  because  of  multiple  scatter 
phenomena,  which  though  linear,  destroy  the  reciprocal  character  of  conventional  solutions  (for 
ideal  media).  Moreover,  the  requirement  of  effective  modeling  inposes  further  complexities:  the 
field  nature  of  the  solutions  must  include  the  role  of  the  transmitting  and  receiving  arrays  whereby 
the  system  sensors  and  signal  processing  itself  are  suitably  coupled  to  the  medium  in  question,  cf. 
Fig.  1.1  ff.  In  addition,  doppler  and  gradient  effects  may  also  be  important  For  short  ranges 
velocity  gradients  in  the  medium  can  often  be  neglected  (i.e.,  Vc  =  0),  but  for  longer  ranges  they 
must  be  carefully  taken  into  account  (Vc  *  0).  Furthermore,  geometry  plays  a  critical  role: 
platform  motion,  beam  size,  structure,  and  direction,  as  well  as  location  in  the  medium,  all  dictate 
potential  interactions  with  the  wave-surface  and  bottom  boundaries,  as  Fig.  2.2  indicates,  along 
with  the  inherently  random  inhomogeneous  effects  of  the  volume  itself  [3],  [4]. 

We  distinguish  two  main  areas  of  interest:  A.  deterministic  inhomogeneous  media ,  where 
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random  variations  are  largely  ignorable  vis-^-vis  the  large-scale  deterministic  variations,  and  B. 
random  inhomogeneous  media,  for  which  significant  media  variations  are  stochastic,  with  and 
without  deterministic  components.  The  former  class  (A)  appears  to  be  so  far  the  one  which  has 
received  the  greatest  attention  with  respect  to  numerical  methods  [7],  [8],  and  [9],  [10],  with 
attention  largely  centered  on  the  purely  spatially  dependent  velocity  gradient  Vc(r,z)  =  c0  +  t(r,z). 
Here,  however,  in  advance  of  current  efforts,  we  shall  direct  our  attention  to  the  more  complex 
situation  (B),  where  in  addition  to  deterministic  variations  of  Vc,  there  are  important  random  ones 
as  well.  In  fact,  our  treatment  will  proceed  in  a  general  canonical  way,  allowing  for  both  the 
deterministic  and  random  variations  in  the  medium,  including  its  important  boundaries. 

1.1  Aims  of  the  Paper 

Before  numerical  procedures  can  be  applied,  however,  it  is  necessary  to  develop  the 
appropriate  analytical  forms  of  solution  from  which  in  most  cases  numerical  results  must  be 
extracted.  The  aims  of  this  paper,  therefore,  are  to  outline  an  analytical-operational  framework  for 
initiating  computational  solutions  to  underwater  acoustic  scattering  problems,  in  particular,  which 
include  the  relevant  ocean  surface  and  bottom  boundaries,  as  well  as  volume  effects. 

We  emphasize  that  this  paper  does  not  intend  to  provide  explicit  computational  procedures 
and  solutions:  it  is  rather  aimed  at  providing  the  "macro-algorithms"  from  and  to  which  appropriate 
software  may  be  designed  and  applied.  Our  aims  include  providing  formal  operational  structures, 
with  some  special  cases  in  which  analytic  solutions  are  obtainable,  and  which  consequently  require 
a  lower  order  of  computational  effort  and  sophistication.  A  collateral  aim  of  this  work  is  to  entice 
interest  in  pursuing  a  computational  formulation,  based  on  the  various  formal  operational 
"solutions"  derived  here  for  the  composite  deterministic  and  random  fields  a(R,t),  namely  <a>, 
cf.  (2.33),  and  higher-order  moments  <aj<X2>,  etc.  Our  approach  here  is  a  "top-down"  approach, 
in  which  we  start  with  as  general  a  formulation  as  possible,  consistent  with  the  broad  class  of 
problems  embodied  in  the  propagation  and  processing  schema  of  Fig.  1.1  and  the  geometries  of 
Fig.  2.1  ff.  We  proceed  then  to  more  specialized  cases.  This  reveals  the  embedding  of  the  latter  in 
the  generality  of  the  former,  and  thus  provides  a  "perspective"  for  specific  medium  models. 

12  New  Elements 

Among  the  new  features  of  die  present  paper  are 

(1)  Formal  operational  formulations  which  can  provide  the  basis  for  a  large-scale 
computational  attack  on  these  four-dimensional  scattering  problems,  Sec.  2. 

(2)  A  feedback  operational  method  for  obtaining  the  (first-order)  equivalent  medium 
renormalization  operator,  cf.  (2.34),  by  means  of  a  converging  series  of  estimates, 
Sec.  3.3. 
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Fig.  1.1  Operational  schematic  of  an  inhomogeneous  (linear)medium,  with  source  and  receiver 
coupling.  Here  Q  is  the  inhomogeneity  operator  associated  with  this  general 
inhomogeneous  (linear)  medium. 


(3)  Some  fully  analytic  solutions  in  important  special  cases,  e.g.,  infinite  or  "weakly- 
bounded"  media,  where  in  effect  there  is  negligible  interaction  between  sources, 
boundaries,  and  among  the  different  scattering  modes  (Surface,  Volume,  Bottom), 
Sec.  2.7. 

(4)  Various  explicit  physical  scattering  models  for  the  inhomogeneity  operator 
including  the  important  wind- wave  scattering  cases,  Eqs.  (2.31). 

(5)  A  hierarchy  of  approximations  for  the  equivalent  deterministic  inhomogeneity  operator 

which  apply  in  different  physical  situations,  Sec.  2.5. 

13  Organization 

Accordingly,  the  present  paper  is  organized  in  the  following  general  but  concise  fashion: 
Section  2  introduces  formal  operational  solutions,  which  form  the  basis  for  the  treatment  of  both 
deterministic  and  random  media.  In  Section  3  various  diagram  methods  and  solution  equivalents 
are  presented,  first  for  the  deterministic  cases  and  then  for  the  many  situations  involving  random 
media  components.  Useful  approximations  are  next  briefly  described  in  this  section,  along  with  a 
possible  method  for  estimating  the  renormalization  operators.  Section  4  completes  the  paper  with  a 
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discussion  of  the  interfacing  problems  between  the  analytic/operational  and  computational 
formulations  by  reviewing  via  the  text  the  main  steps  leading  up  to  the  as  yet  unachieved  software 
realizations. 

2.  Formal  Operational  Solutions 

Accordingly,  the  above  discussion  indicates  the  critical  need,  at  some  point,  for  numerical 
procedures  to  obtain  the  needed  quantitative  description  of  the  received  acoustic  field,  for  whatever 
operational  or  oceanographic  purpose  at  hand. 

What  arc  the  needed  descriptors  of  this  received  field?  This  will,  of  course,  depend  on 
whether  or  not  there  arc  significant  random,  as  well  as  deterministic  inhomogeneities  in  the 
medium.  Here  we  shall  assume  a  (linear)  combination  of  both  types,  which  permits  us  to  focus  on 
one  or  the  other  as  limiting  cases.  Thus,  in  general  we  can  write  formally  for  the  inhomogeneity 
operator  0.  cf.  Fig.  1.1 

Q=$(h)  +  $(I),  (2.1) 

where  0®  are  respectively  associated  with  the  deterministic  and  random  parts.  For  example, 
Qw  usualy  dominates  in  a  nonzero-gradient  ocean  with  small  surface  interactions  (cf.  some  of  the 
current  treatments  [7]— [10]),  while  Q®  is  most  significant  in  surface  or  bottom  scatter  situations,  at 
comparatively  short  ranges  and  high  frequencies,  where  Vc  =  0  often,  cf.  [1],  [4]. 

Now,  since  the  medium,  including  its  boundaries,  is  both  spatially  and  temporally 
stochastic  in  general,  "solutions"  to  the  appropriate  propagation  equations  for  the  input,  incident, 
and  scattered  fields  are  necessarily  statistics  of  these  fields,  namely,  means,  variances, 
covariances,  higher-order  moments,  and  finally,  various  orders  of  the  probability  distributions  of 
these  fields.  Thus,  if  a(R,t)  represents  the  (acoustic)  field  at  point  (R,t)  in  this  random, 
inhomogeneous  medium,  one  seeks  «x>,  «Xi<X2>  (=  <a(R1,t1)a(R2,t2)>),  <a1a2a3>,  etc.,  and 
more  completely,  w^ctj),  w2(a1,02),  etc.  The  lower  order  moments,  for  example,  «X!>, 
<a!0t2>,  are  needed  in  obtaining  the  scattering  cross-sections  of  the  "illuminated"  wave-surface  or 
bottom,  cf.  Fig.  2.1,  while  w^c^)  is  required  for  signal  processing  (cf.  Part  II  of  [4]).  What 
makes  the  latter  critically  non-trivial  is  the  fact  that  a  is  often  strongly  nongaussian  (cf.  VI  of  [4]).. 
Even  in  the  cases  where  the  Central  Limit  Theorem  applies  (cf.  Sec.  7.7  of  [5],  for  example),  so 
that  the  field  is  gaussian,  one  still  needs  the  lower-order  moments  <cx>,  <aja2>,  in  order  to 
specify  the  gaussian  process,  a(R,t). 

2.1  Formal  Development 

Clearly,  for  realistic  solutions,  many  complicating  effects  must  be  included,  which  almost 
always  at  some  point  makes  a  computational  effort  essential.  Initially,  one  starts  with  a 
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Fig.  2.1  Schematic  scatter  channel  geometry  of  source  (T),  medium  (V),  and  boundaries  (SJ1), 
with  receiver  (R),  D  =  unscattered,  or  "direct"  path,  and  n  is  the  order  of  the  interactions. 
Sufficiently  high  frequencies  are  assumed  to  permit  illustrative  ray  paths. 

propagation  equation  which  at  least  implicitly  contains  the  major  physical  phenomena  governing 
the  generation  and  development  of  both  the  scattered  and  unscattered  fields.  Formally,  for  the 
linear  media  assumed  here,  the  propagated  field  a(R,t)  obeys  a  partial  integro-differential  equation 
of  the  form 


Here  Gy  is  the  source  function,  e.g.,  GT  =  6xSin,  cf.  Fig.  1.1,  and  is  a  linear  (scalar)  partial 
differential  operator,  where  is  associated  with  the  homogeneous  portion  of  the  medium.  In 
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(2.2),  0  is,  generally,  a  (scalar,  linear)  integro-differential  scattering  operator,  which  describes  the 
interaction  of  the  incident,  i.e.,  homogeneous  field  with  the  differential,  or  local,  scattering 
elements  of  the  inhomogeneous  portion  of  the  medium.  (The  possible  integral  operator  component 
of  0  stems  from  local  scatterer  doppler  effects,  for  example.)  Here,  b.c.'s  =  boundary  conditions, 
and  i.c.'s  =  initial  conditions.  Since,  formally  from  (2.2)  a  =  (£,(°>  -  Q)*1(-Gr)lb.c.+i.c.,  from  a  = 
T<JP(-Gt).  cf.  Fig.  1.1,  one  has  directly 


T^  =  (t(°>-Q)-1  (2.3) 

for  the  medium  operator  T^,  now  expressed  in  terms  of  the  homogeneous  and  inhomogeneous 
components  of  the  medium. 

Because  of  the  usual  statistical  nature  of  Q,  Eq-  (2.2)  is  really  an  ensemble  of  equations, 
which  is  called  the  (associated)  Langevin  Equation  [cf.  Chapter  10  of  [5]),  which  can  be  written  as 


{(£(0)  -  0)a  =  -  GT  +  (b.c.'s  +  i.c.'s)). 


(2.4) 


where  the  brackets  {  )  emphasize  the  ensemble  character  of  this  Langevin  equation.  In  subsequent 
discussion  the  brackets  may  be  omitted  for  convenience,  it  being  noted  that  Q,  and  therefore  a,  are 
stochastic.  (See  Sec.  IIB  of  [4]  for  remarks.) 

An  important  example  for  acoustic  propagation  in  underwater  media  is  the  case  of 
inhomogeneous  absorptive  media  [3],  for  which  (2.2)  is  specifically  [11] 


where 


{(■ 


+  *x(r,t) 


atj 


2  Id  +  e(R,t)] 


an 

dt*  J  a 


-Gj 


(2.5a) 


£<°)  =  fi 


+  tox 


iWii 

dt)  eg  dt2 


a  e  d 

(2-5b) 

with  Tx  =  x0x  (1  +  Yx[*  ,t])  generally.  When  xox  =  0,  (2.5a)  reduces  to  the  more  familiar  "extended" 
Helmholtz  equation 


{ 


1 


^11 


V2_-L[(i+e(R,t)]—  I  <x  =  -GT 

4  3t2  J 


(2.6) 
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with 


where 


A-i.  <fL. 


t(°)  B  V2  - 


jl  # 

cj  dt2’ 


(2.6a) 


R  =  ixx  +  iyy  +  izz , 


(2.7) 


and  rectangular  coddinates  are  assumed  throughout,  so  that  V  •  V  =  V2,  etc.  Specifically,  co  = 
(constant)  wavefront  ?ed,  e  embodies  the  effects  of  velocity  gradients,  internal  wave 
phenomena,  and/or  (v  ok)  local  turbulence,  while  xx  represents  the  effects  of  relaxation  absoprdon 
(due  to  Mg2S04  and  other  salts). 

Similarly,  in  the  more  restricted  situations  where  (2.6)  is  replaced  by  a  purely  deterministic 
model  in  which  e(R,t)  — >  e(R)  is  nonrandom,  this  Helmholz  equation  can  be  approximated  (at  a 
given  frequency  f0)  by  a  basic  form  of  the  parabolic  equation  (PE),  which  is  espressed  by  (Eq. 
(2.34),  [8]),  in  cylindrical  coordinates  [R  =  (r,z)  =  (r,0,z)]: 


^=(A  +  BS“°;  “0=“o(R>> 

where  specifically 


(2.8) 


A  =  |  ko  [n2(r,z)  -  1];  Bs^;  ko  =  <Uo/c0  =  2jtfo/c0 ,  (2.8a) 

n2  «  (1  +  e)  .  (2.8b) 

Comparing  (2.2)  and  (2.8)  gives  at  once  here 

=  0  =  Ael  =  [-ikoe/2]t  (2.9) 

dzl  dr 

where  1  =  unit  operator,  1c  =  C,  etc.,  and  (Xq  =  0  is  an  equation  of  the  Schrodinger  type  since 
B  is  purely  imaginary. 

The  formal  solution  of  the  Langevin  equation  (2.4)  is  obtained  from  the  equivalent 
operational  form 


(f  -  )ClQ)a  =  lft(-  Gx) ;  a  =  aH  +  04;  aH  =  1&(-  GT), 

8 


(2.10) 
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where  aH  is  the  homogeneous  (i.e.,  non-scattered)  component  of  the  total  field  a,  with  aj  the 
associated  inhomogeneous  (scattered)  field.  A  is  £,(°H,  with  the  boundary  and  initial  conditions; 
in  fact,  A  is  the  integral  Green's  function  operator 

— 

A(R,t  I  R',f)  =  -  Jdf  JdR'g(R,tlR',t’)  ( )R  t-  (2.11) 

—  (-) 

obtained  from  f/°)g  =  -5(R-R')  8(t-t').  The  formal  operational  solution  of  (2.10)  now  becomes 

(2.12) 

Herer|  is  the  field  renormalization  operator  (FRO),  and  (2.12)  represents  the  Feedback  Operational 
Solution  (FOS)  for  a,  c.f.  remarks  following  Eq.  (2.35)  et  seq.  Expansion  of  (2.12)  yields  the 
Perturbation  Theoretical  Series  Solution  (PTSS) 


(2.13) 

where  =  [■nq...'n]n  is  the  n'Mteration  of  rj  and  where  the  operator  formalism  used  here  is 
represented  explicitly  by 

(A)(n)  =  A(»>  £  A(R,t  I  R',t')  A(R',t'l  R",t")  .  .  .  A(R(^D,t(n-l)  |  R(n)t(r))  (2.13a) 

which  is  generally  distinct  from  the  direct  operator  product 

A"  £  A(R,t  I  R’,0  A(R,t  I  R",t") . . .  A(R,t  I  R("),t(n>)  (2.13b) 

and  its  more  general  form 

Ai . . .  A"  =  A(R,t  I  R',t') .  .  .  A(Rn,tn  I  R(n),t<n)).  (2.13c) 

We  have  also  the  conventions 
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Ml  =  A(Ri,tilR,,t’)^(Ri,tilR",t") 
M2  =  A(Ri,tilR,,f)^(R2,t2IR",t") 


}  (2.13d) 


so  that  clearly  A1&2  ^  &2&1.  i.c.,  Ai,  62  do  not  commute,  unless  the  operand  is  symmetric. 

22  Some  Homogeneous  Methods  and  Results 

For  the  usual  cases  where  there  is  negligible  interaction  of  the  medium  and  boundaries 
upon  the  source,  we  can  replace  A  by  &»,  where  is  the  integral  Green’s  function  operator  for 
the  associated  infinite  homogeneous  medium.  Then,  also,  we  have  t\  — »  cf.  (2.12). 

Next,  to  proceed  further  it  is  clear  from  (2.12)  that  we  must  determine  the  homogeneous 
field  cch,  =  0,  or  0  =  0.  now.  In  the  light  of  the  fact  that  media  boundaries  or  interfaces  can  be 
equivalently  considered  as  (distributed)  inhomogeneities  in  an  otherwise  homogeneous  medium, 
cf.  Sec.  2.4  ff.,  to  obtain  solutions  in  the  homogeneous  cases  we  must  therefore  consider 
unbounded  media,  along  with  the  condition  Q  =  0.  More  precisely,  by  "unbounded"  we  mean  the 
infinite  spatial  domain  exclusive  of  sources. 

The  general  method  of  solution  here  (including  bounded  media)  involves  two  steps: 


I.  Obtain  the  associated  Green's  function,  g,  from  (2.2),  (0  =  0); 

II.  Construct  the  corresponding  GHP  (Generalized  Huygen’s  Principle)  [20], 


and  thereby  obtain  the  field,  a.  For  unbounded  homogeneous  media  this  is  easily  done.  Thus,  the 
Green's  function  g«,  is  found  from  the  defining  relation 


t(0)gao  =  -5(t-t')8(R-R'), 


(2.14a) 


g(R,t  I  R',t')oo  =  -  iftwSt't-Sinr , 


(2.14b) 


where  Kfcoo  is  the  inverse  (integral)  operator  to  £(°)  +  i.c.’s,  with  =  £P).  Integrating  over 
(R',t')  shows  at  once  that  iCt*.  and  g«  are  related  here  by 
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where,  by  causality,  g*  =  0,  t  £  t'-,  which  represents  the  initial  condition  (i.c.)  on  the  "point"- 
source  in  (2. 14a).  Thus,  is  now  the  Green’s  function  operator ,  with  the  kernel  goo. 

The  desired  field  a(R,t)oe  is  found  directly  from  (2.2)  (Q,  b.c.'s  =  0),  viz.: 


(2.16a) 


(2.16b) 


(2.16c) 

where  Gj  =  o,  t*  <•  to-.  R  *  Vj:  the  source  function  is  initiated  at  t'  =  to  and  occupies  the  domain 
Vt(<°°).  The  relations  (2. 16b,c)  are  the  simplest  form  of  GHP  here  for  homogeneous  unbounded 
media  (outside  the  source  domain). 

Various  transform  equivalents  of  prove  useful.  These  may  be  expressed  as  follows, 
from  (2.16b): 

•o 

=  {  fdt*  |Y0(R,s ! R't')«, ( )R',t' dR'}  h-Js-1  {Y0,oo}  (2.17a) 

—  (-) 


— {  Jdt'  jY0(k,s  I  R’t')oo  ( )R',t'  dR' }  s  /  0,«)  (2.17b) 

-oo  (oo) 


£(0)aoo  =  -  Gt,  or 


—  Jdt  6x(gooSin)  —  jTATSindt , 


=  -^kJs-1  {ms)oo  Jdt’  Je_ik,R_st'  ( )R>1'  dR'}  (2.17c) 

— °  (~) 

where  (2.17a-c)  define  the  transformed  operators  £0,oo,  5o,~  (=  $o(k,s  I  R',t')«,)  respectively, 
and  where  {A»}  =-^0.00.  ^R_1  =  -  3o(k,s  I  R',t')  are  the  corresponding  inverses. 

The  general  transform  relations  are  explicitly 


a(t-t')  =  Jr1  {e-st’A(s)}  = 


(2.18a) 


where  s  is  a  complex  (angular)  frequency,  and  provided  all  the  singularities  of  A(s)  lie  in  Re(s)  < 
d;  a  =  0,  t  <,  t'-  is,  of  course,  an  initial  condition  (i.c.).  The  inverse  of  (2.18a)  is 
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c_st  A(s)  =  y\  { a(t— t*) }  =  Je-sta(t-t')dt  =  e-sl'  Je-*a(t)dT,  Re(s)>0.  (2.18b) 

In  the  steady-state  regime  (t'  — » -»)  we  choose  d  =  0  and  use  s  =  ico  =  2jtif,  where  (co,f)  are  real 
frequencies,  with  an  appropriate  indentation  of  the  contour  (Bri,  d=0)  for  any  singularities  on 
Re(s)  =  0.  The  contour  BRi  is  a  Bromwich  contour  [26].  Similarly,  Jr-1,  7k  are  spatial 
Fourier  transforms,  e.g., 


B(k)  =  ^-Mb(R)}=  Je-*-Rb(R)dR; 

[R1 


b(R)  =  7k  (B(k))  = 


]eik*RB(k) 

(k) 


dk 

(2rc)3 


(2.18c) 


where,  as  before,  cf.  (2.7),  k  =  f*kx  +  tyky  +  tzkz  (=  2rcv)  is  a  vector  wave  number  (and  v  is  a 
vector  spatial  frequency).  Multiple  transforms  of  quantities  like  C(R,t)  are  handled  in  the  same 
way,  cf.  (2.17b,c). 

For  those  unbounded  media — unbounded  in  the  sense  that  a  propagating  wave  front  never 
encounters  a  bounding  surface — the  evaluation  of  the  Green's  function  goo  from  (2.14a)  and  the 
general  field  from  (2.16)  is  readily  accomplished.  From  (2.14a)  one  has 


Jr-1  £(°)  { 7s ~l7k  {go.-} }  =  (2.19a) 


/.  gof«  =  go(k,s  I  R',t')oo  =  -Lo(k,s)-1e-ik,R,-st' ,  (2.19b) 

where  Lo(k,s)  is  the  result  of  applying  tP)  to  Js-1^k>  remembering  that  Ve**11  =  ikeik'R, 

V^ik’R  -  -k*keik*R  =  -k2eik*R,  —  est  =  sest,  etc.  In  fact,  one  has 

dt 


-Lptk.s)-1  =  %(k,s)oo 


(2.20) 


cf.  (2.17b,c)  so  that,  from  (2.19b), 

g(R,t  I  RV)oo=  Js"1  Jc  (go.-)  =  Je^-O—  Jeik*(R-R'){-L0(k,s)-1}  .  (2.21a) 

Bri  2rti  [k)  (2rt)3 

This  can  be  put  in  the  more  convenient  equivalent  form  for  evaluation  with  the  help  of  spherical 
coordinates:  k  =  (k,  w  [=  cos  0],  0),  viz.: 
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g(R,t  i  R't'). 


,=  Jes(t_t')  — 
Br!  2ia 


,  l  2k 

-  "  ^  /k2dk  Jd w  •  J[-L0(k,s)_1]eiklR_R’,w  d<(> 


(2.21b) 


A.  Conditions: 

Various  conditions  and  properties  of  (2.21b)  are  now  noted: 


(i).  "Spatial  causality this  is  the  radiation  condition  [15] 


( 


lim 

R— »oo 


V 


dj 

9R 
(.  Rik  . 


>  + 


R 

1  a' 

\ 

r 

• 

c°  dt 

► 

1 

* 

R 

s_ 

^  c0  ' 

LYq.o* 

=  0,  R  =  IRI , 


J 


(2.22) 


which  ensures  that  only  outgoing  waves  from  the  source  (eVt)  are  propagated.  This  spatial 
causality  condition  is  used  in  the  specific  evaluation  of  go>0o,  (2.21b),  by  appropriate  modifications 
Ck  of  the  k-contour  (-“,«>)  to  include  only  singularities  in  the  complex  k-plane  which  yield  the 
required  outgoing  waves. 

(ii).  "Temporal  causality this  is  the  well-known  condition  that  propagation  cannot  occur 
before  the  source  initiating  it  is  activated,  cf.  (2.15),  (2.18),  or  more  precisely,  that  events  cannot 
occur  (macroscopically)  at  some  point  P(R,t)  before  the  source  has  been  initiated  and  the  field 
generated  in  the  medium  has  reached  the  point  in  question.  For  causality  in  time  one  has  the  Paley- 
Wiener  criterion  [5],  pp.  96, 97, 102  and  Refs.  Here  (s  — >  io>): 


rioglY0(Rja)IR%  P 

J  1+0)2 

o 


dco  <  oo  , 


(2.23) 


where  specifically 


Y0(R,slR')oo  =  { g(R,slR',t')ooest'}  =  Y0(R,slR',t’)«est’,  (2.23a) 


,  1  2k 

/k2dk  Jdw  j[-L0(k,s)]-l]e*IR-R^d<)>, 


(2.23b) 


from  (2.21b),  (2.17a).  This  condition  is  obeyed  generally,  from  the  initial  condition  on  g,  cf. 
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(2.15).  The  combination  of  (222)  and  (223)  is  the  full,  space-time  causality  condition,  which 
when  axiomatically  invoked  (as  it  is  here)  ensures  that  only  time-related  solutions  of  (2.16)  are 
possible. 

(iii) .  Reciprocity :  the  point  source  [in  (2.14a)]  and  the  observer  [at  P(R,t)]  are 
interchangeable  in  the  Green's  function,  provided  temporal  causality  [(ii)  above]  is  also  obeyed. 
Thus,  one  has  R'  -» R,  R  -»  R'  and  t  -» -t',  t'  — » -t,  so  that  g(R',-t'IR,-t).  From  (2.21)  it  is  at 
once  evident  that  reciprocity  for  these  homogeneous  media  is  obeyed. 

(iv) .  "Regularity  at  infinity [15]  this  ensures  that  igeol  (£  <»),  i.e.,  is  bounded  as  R  — » »; 
(this  is  not  necessarily  guaranteed  by  (2.22)).  This  condition  is  a  reflection  of  the  fact  that  the 
(point-)  source  has  bounded  energy  and  consequently  its  field  (goo)  must  have  also,  all  R. 

Not  too  surprisingly  from  the  above,  we  may  now  regard  these  homogeneous  linear  media 
as  being,  in  effect,  linear  space-time  filters.  In  the  language  of  modem  circuit  theory  the  Green’s 
function,  goo,  represents  the  medium's  weighting  function',  Y0too,  %»  correspond  to  the 
frequency-,  and  frequency-wave  number  system  functions,  or  "filter"  spectral  densities,  of  this 
causal,  space-time  filter,  where  the  governing  space-time  causality  condition  is  the  combination  of 
(i)  and  (ii)  above,  cf.  (2.22),  (2.23).  However,  unlike  conventional,  "lumped-element"  circuit 
models,  these  "circuits"  (i.e.,  the  medium)  are  in  effect  composed  of  a  continuum  of  elements, 
distributed  throughout  space.  The  dynamical  ("circuit")  equations  of  the  latter  are  partial 
differential  equations,  obeying  boundary  conditions  (here  as  R  -* «»,  cf.  (2.22))  as  well  as  initial 
conditions,  while  the  former  are  limited  to  ordinary  differential  equations,  subject  only  to  initial 
conditions;  (both  obey  appropriate  stability  conditions,  (iv),  for  example),  [27],  Further  critical 
differences  appear  when  the  medium  becomes  inhomogeneous,  and  random  [cf.  Sec.  2.4  ff.]. 

For  the  Helmholtz  case  (2.6),  (2.6a),  (^  =  0,  it  is  easily  seen  that  (2.1 1)  becomes  (2.15), 
and  that  following  the  methods  described  in  (2.16)— (2.18)  one  gets  specifically  for  with  the 
Green's  function  now 


gDc  =  5(At-p/c0)/47tp  ;  p  =  IR  -  R'l,  At  =  t-t'.  (2.24) 

Similarly,  it  is  also  found  [3]  that  when  0  =  0,  A  ->  Aoo,  the  solution  of  (2.5a)  reduces  to 
g(R,t  I  R’,t')oo  =  FkFs-l  {  ft,.-)  =  Fs-l  {?  o,~),  (2.25a) 

where  F,  etc.  denote  Fourier  transforms,  with 

ft>,~  =  Jo(k,s)  f dt'  Je-ik*R-sl'(  )R- t.dR’  =  - Fr-1  Ft  {id  »),  (2.25b) 

J  M 
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%(k,s)«  =  [k2(l+toXs)  +  S2/CC2]-1 . 


(2.25c) 


It  can  be  shown  that  here  [6] 


<>  _ 

A  O.oo  — 


e-(p«/CoXlvioX*)-t/2-*f 

47Cp(l+toXs) 


(2.26a) 


and  that  now,  explicitly, 


goo  =  -1—  (  1  -  [  r— P1  c*A^t°xiFi(-l-2n;j  -n;At/toX)  ) 

47tptox  l  (2n+l).  2  J 

(2.26b) 

for  the  Green's  function  of  the  associated  unbounded  medium,  where  iFj  is  a  confluent 
hypergeometric  function. 

Finally,  in  the  case  of  the  PE  approximation  (2.8)  above,  one  can  show  that  the  Green’s 
function  associated  with  (2.9)  is 


g(r,zlr',z')oo  =  e^2^,  Ar  =  r-r'. 

2rcAr 

With  a  gaussian  beam  for  the  source  (along  z)  for  these  narrow  angles,  e.g., 

Gt  =  Aoe-z2/a2 , 


(2.27) 


(2.28) 


the  associated  homogeneous  field  Oqh  is  found  to  be  [19] 


Ooh(At,z)  =  Jg(r,zlr',z')0O  GtU'Mz’  , 


(2.29a) 


akpAp  T  k2a2z2  .  f  2ArkoZ2  .  2Ar  VI 

~  (4Ar2+k2a4)1^2  CXP  [”  4Ar2+k2a2  +  *  (4Ar2+k2a4  "  Un’  k^Jj  ’  (2,2%) 

which  shows  the  expected  (Ar)-1  fall-off  as  Ar  — » 

23  The  FOR 

The  result  (2.12)  above  can  be  expressed  formally  in  terms  of  a  generalized  feedback  loop, 
as  shown  in  Fig.  2.2.  Thus,  (2.12)  can  be  regarded  as  an  extension  of  conventional  control 
theory,  now  to  deterministic  inhomogeneous  media  for  the  moment,  cf.  Sec.  2.4  ff.  We  call  Fig. 
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2.2  the  feedback  operational  representation  (FOR)  for  the  FOS,  (2.12).  Equation  (2.13)  above  is 
the  associated  PTSS. 


e(R.t) 


Fig.  22  Feedback  Operational  Representation  (FOR)  for  the  propagated  field  a,  Eq.  (2.12  ),for 
linear  inhomogeneous  media. 


2.4  The  Inhomogeneous  Operators  Q 

These  operators  embody  the  various  inhomogeneities  characteristic  of  the  medium, 
including  its  boundaries.  Thus,  for  a  typical  propagation  scenario,  sketched  in  Fig.  2.1,  we  may 
express  Q,  (2. 1 ),  in  more  detail  for  these  linear  media  as  a  combination  of  scattering  modes: 

Q =  Os  +  Qv  +  Qb  =  [$?  +  0?]  +  [ft?  +  +  [#£>  +  0^  (2.30) 

where  S,  V,  B,  respectively,  denote  the  surface,  volume,  and  bottom  in  this  instance,  and  (h),  (I), 
designate  the  deterministic  and  (purely)  random  components  of  each  associated  operator,  cf.  (2.1). 

As  examples,  we  have  explicitly  from  the  various  acoustic  propagation  models  indicated  in 
Section  2.1,  and  for  interfaces  cf.  [1],  Eq.  (15);  [18]: 

•(it-1r)Vr|£sj=R0S  h»(iT-fR)s/co;  s  =  ico,  (2.31a) 

for  the  Kirchoff  approximation  to  the  surface  and  bottom  operator,  with  £  =  wave  or  bottom 
surface,  cf.  Fig.  2. 1 .  Here  R0,S  are  the  plane  wave  reflection  and  shadowing  functions,  and  A  is 
the  (inwardly  drawn)  normal  to  £.  Also,  Ij,  Tr  are  unit  vectors  in  the  directions  of  the  incident 
and  scattered  wavefronts.  A  far-field  or  Fraunhdfer  condition  is  assumed.  For  volume 
inhomogeneities  we  have 

Qv  =  Eq.  (2.5b),  or  Eq.  (2.6a);  Eq.  (2.9)  for  the  PE  approx.  (2.8),  etc.  (2.31b) 
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The  role  of  T|(n)  in  (2. 1 3)  is  to  generate  multiple  scatter  components: 

ft2)  =  (fooQ)(2>  =  Life  +  Qv  +  ^b)2} 

=  £{$?  +  QsQv  +  ftvfts  +  •••)  (2-32) 

shows  the  various  second-order  multiple  scatter  contributions  among  the  various  types  of 
inhomogeneity.*  In  most  underwater  propagation  problems  inter-modal  scattering  can  be 
neglected,  e.g.,  terms  in  QsQv.  etc.  are  ignorable  vis-k-vis  the  intra- modal  contributions  Q(2), 

*  which  condition  can  considerably  simplify  the  subsequent  calculations.  In  any  case,  the 
structure  (2.30),  or  "M-form,"  appearing  in  the  FRO,  T),  provides  a  useful  "bookkeeping"  of  the 
various  possible  interactions. 

With  inhomogeneous  (Q  *■  0)  media,  exact  conventional  techniques  fail,  largely  because  the 
medium  is  not  reciprocal  (cf.  (iii),  Sec.  2.2]  and  there  is  no  general  way  of  applying  and  evaluating 
the  initial  conditions  over  the  various  volume  integrals  which  appear  in  the  development  of  the 
Generalized  Huygens  Principle  (GHP)  [20]  when  Q  *  0.  [For  examples  of  the  GHP  see  Sec.  7.3, 
[20],  and  pp.  512,  513,  [28].]  Even  in  the  simplest  cases,  which  are  bounded  homogeneous 
media — the  boundary  (X)  being  the  inhomogeneity  here — the  usefulness  of  the  results  is  severely 
limited  by  the  difficulty  in  obtaining  the  corresponding  Green's  functions,  gx,  except  for 
boundaries  compatible  with  separable  coordinate  systems  (for  which  complete  sets  of 
eigenfunctions  can  be  obtained).  A  second  difficulty  is  that  the  resulting  GHP  is  basically  an 
integral  equation  in  the  desired  scattered  field  (i.e.,  scattered  by  the  boundary  surface  X),  a,  unless 
it  makes  physical  sense  to  specify  the  field  h.  priori  on  X,  an  unlikely  situation  in  most  applications. 
A  third  difficulty  is  simply  the  fact  that  these  simple  cases  are  not  the  ones  of  principal  interest 
Finally,  a  fourth  complication  is  that  0  is  often  an  integral  operator,  e.g., 

&R,t  I  R’,t')  =  Jdt’  jQ(R,t  I  R’,f)  (  )R',f  dR'  ,  (2.32a) 

cf.  (2.15)  for  ifao*,  for  example,  formally,  expressing  the  physical  fact  that  the  effects  of  the 
inhomogeneities  may  not  be  local  in  their  interactions. 

Another  important  reason  why  conventional  techniques  fail  in  the  treatment  of 
inhomogeneous  media  is  the  fact  that  these  media  are  space-time  variable  in  their  inhomogeneities, 
e.g.,Q  =  Q(R,tl..  .),  cf.  (2.5),  (2.6).  Such  media  do  not  generally  support  space-  and/or  time- 
harmonic  solutions,  exceptions  being  Q(RI ...),  which  permits  time-harmonic  solutions,  and 

*  Of  course,  one  has  an  appropriate  preceding  each  $ 
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frtl...)  ,  which  allows  space-harmonic  solutions.  Even  the  standard  perturbational  and  variational 
techniques  of  "classical"  approaches  (Chapter  9,  [20])  are  only  applicable  in  the  case  Q(RI  ...), 
when  @  is  deterministic,  which,  of  course,  is  often  not  the  physical  problem  in  ocean  and 
atmospheric  media,  whose  local  and  distributed  properties  are  explicitly  space-  and  time- 
dependent.  Moreover,  even  these  methods  break  down  when  the  medium  is  random  (i.e., 

&RI  -  .)  is  stochastic).  In  fact,  all  the  usual  deterministic  methods  become  nugatory  when  the 
medium,  and  its  boundaries,  are  random  and  moving. 

Even  though  exact  analytical  and  standard  (i.e.,  nonrandom)  perturbational  techniques  fail, 
various  useful  approximate  methods  are  available.  The  best-known  of  these  is  ray-tracing  (a 
comparatively  high-frequency  technique  [21],  [22];  [29],  Parts  HI,  IV),  and  modal  analysis  (a 
comparatively  low-frequency  one  [21]).  An  important  recent  extension  of  these  ideas,  combined 
with  path-integral  methods  and  super-eikonal  procedures,  has  been  developed  and  applied  to  the 
study  of  sound  speed  fluctuations  in  nonisotropic  oceans  by  Flattl  et  al.  [22].  All  these  methods 
are,  of  course,  necessarily  approximate  and,  in  fact,  skillful  approximation  is  the  successful  heart 
of  these  approaches.  In  our  present  approach  to  scattering  problems  we  seek  to  preserve  the  exact 
formal  solutions  in  a  canonical  development,  reserving  approximations  to  the  last  possible 
moment  This  has  the  conceptual  advantage  of  explicitly  identifying  the  disposable  portions  of  the 
general  result,  in  the  contest  of  the  particular  problem  ultimately  to  be  studied.  Approximations  are 
inevitable:  we  attempt  to  postpone  them  as  long  as  possible.  Accordingly,  in  this  spirit,  we  seek 
formulations  which  may  lead  to  manageable  computational  (and  occasionally  analytic)  solutions  of 
(2.2),  (2.4),  as  well  as  using  older  methods. 

25  Statistical  Solutions:  Dyson  and  Bethe-Salpeter  Equations 

In  view  of  the  statistical  character  of  the  scattering  mechanisms  in  most  cases  of  physical 
interest,  as  already  noted  above  at  the  beginning  of  Sec.  2  and  in  Sec.  2.1,  t\,  (2.12),  is  a 
stochastic  operator.  Thus,  the  "solutions"  of  the  Langevin  equation  (2.4),  as  mentioned  above,  are 
the  various  moments  of  the  total  field  a,  viz. 

<a>  =  aH  +  <  X^(n)>aH.  (2.33) 

11=1 

etc.  for  «Xi02>,...,  with  (J&  -*  ($«  now,  cf.  (2.14a).  Thus,  to  get  the  first  moment  of  a  all 
(n£l)-order  moments  of  the  FRO  are  required.  Similarly,  all  (n,n'£l)-moments  of  <q(n>q(n')>  are 
needed  for  <a1a2>,  etc.  This  is  clearly  very  cumbersome  and  also  presents  questions  of 
convergence-in-probability  (CIP),  e.g.,  Ilqll  <  1  — > <llr|ll>, ...  <llq(n)|l>,  ...  <H?). 

A  way  around  these  difficulties  is  provided,  however,  by  introducing  an  equivalent 
deterministic  relation,  obtained  from 
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<Qa>  s  Oi(d)<a> ;  T)i(d)  =  Oi(d)  •  (2.34) 

Rewriting  (2.12)  as  (1  -  q)a  =  aH  and  averaging,  using  (2.34),  then  gives  the  (deterministic) 
Dyson  equation 


<a>  =  (f — =  aH  +  ^  tj  (dXn)aH  t 

n=l 


(2.35) 


where  now  is  the  Equivalent  Deterministic  Field  Renormalization  Operator  (EDFRO). 

The  evaluation  of  the  nonrandom  equivalent  inhomogeneous  operator  Qi<d)  proceeds  from 
the  relation 


<V-W>  «h  =  <^Ar->  aH 

-  Tloo)^  ^(1  -  r\oo)^ 

obtained  from  (2.12)  averaged  and  (2.35).  Expanding  ^(4)  by 

$i<d>  =  £A$(R,t  IR’,0  H  \  !  Q1W(R,t  IR’.t’)  (  )R',t'  dR'dt 

m=0 

and  applying  (2.37a)  to  (2.36),  with  the  requirement  that  each  term  in  the  resulting  PTSS  of  both 
sides  of  (2.36)  be  of  the  same  operator  order,  give 

~  ~  m21.  (2.37b) 

j=i 

Thus,  0i(d)  is  always  an  integral  operator,  cf.  (2.37a).  When  Q  is  purely  a  fluctuation  operator, 
e.g.,  <Q>  =  0,  which  often  occurs  in  volume  scattering  when  boundary  effects  can  be  ignored, 
usually  for  geometrical  reasons,  (2.37b)  simplifies  considerably.  Thus,  for  example,  we  find  that 
then 


(2.36) 


(2.37a) 


=  0;  A<P  =  <§&£>;  A(f>  =  <Q&2>fy2>>;  ...  .  (2.37c) 

The  extension  to  the  second-order  cases  is  made  in  a  similar  way  ([12],  Sec.  4.2).  For  the 
second-order  moment  <ai(X2>,  which  is  part  of  the  general  statistical  solution  to  the  generic 
Langevin  Equation  (2.4),  we  get 
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<«ia2>  =  (f  -nigLr1  d  -tiff  Tiff)  «*i>  <a2>  , 


(2.38) 


where  rjff  =  etc.,  and 


(2.38a) 

and  <ai>  =  <a(ri,ti)>,  etc.  Here  t)  is  the  second-order  EDFRO,  cf.  (2.35)  et  seq.,  and 
may  be  determined  by  an  extension  of  the  approach  for  above  (2.36)  et  seq.  It  is  found  ([12], 
Sec.  4.2)  that  in  the  frequent  cases  where  <Q>  =  0, 


(2.38b) 

Writing,  cf.  (2.37a), 

OO 

=  £A(2)(R1,t1;R2,t2 1  Rr,tr;R2\t2’), 

m=0 

(2.39) 

and  using  this  in  the  analogue  of  (2.36),  comparing  terms  of  similar  order  in  the  resulting  PTSS 
expansions  gives  finally 

m-0: 

A(o)  =  <0iCfe>  -  <Qi>  <$2>  +  ; 

(2.39a) 

m  =  l: 

A(f =-[<Qi£k>-  <Qi>  <&>]  +  [Ai.~<0i>  +  ^2.~<Q2>] ; 

(2.39b) 

m-2: 

=  etc., 

(2.39c) 

In  the  important  cases  where  <Q>  =  0,  e.g.,  0  is  a  purely  fluctuation  operator  as  is  often  the  case 
for  infinite  media,  all  the  odd-order  terms  in  (2.39)  vanish,  cf.  (3.15). 

Equation  (2.38)  is  a  form  of  Bethe-Salpeter  Equation ,  which  is  the  second-order 
generalization  of  the  Dyson  Equation  (2.35),  and  like  it,  is  now  completely  deterministic.  Still 
higher-order  formulations  for  <aia2a3>, ...,  <aia2...am>, ....  may  be  derived  in  the  same 
formal  way.  These,  however,  appear,  in  most  instances,  to  be  too  complex  to  be  useful. 

Finally,  we  note  again  (cf.  Sec.  2.2)  that  an  important  feature  of  the  governing  Langevin 
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equation  (2.12),  (2.13),  and  the  Dyson  equation  (2.35),  is  that  they  also  embody  *  feedback 
representation,  as  shown  in  Fig.  2.2.  (For  the  Dyson  equation  (2.35),  one  replaces  Q  by  Qi(d) 
and  a  by  <a>.)  Thus,  these  scattering  problems  may  be  regarded  as  an  extension  of  modern 
control  theory,  where  now  the  (linear)  ordinary  differential  equations  of  the  control  format  arc 
replaced  by  partial  differential  equations  of  propagation  in  linear  (inhomogeneous)  random  media. 
Similarly,  initial  conditions  of  the  former  are  replaced  by  both  boundary  and  initial  conditions  for 
the  latter.  This  suggests  that  concepts  and  methods  already  available  from  (linear)  control  theory 
can  be  helpful  here  in  achieving  quantitative  solutions  for  <a>,  <a1a2>,  etc. 

2.6  Related  Approaches:  The  Transport  Equation 

Instead  of  working  with  the  Dyson  equation  (2.35),  we  can  achieve  formal  solutions  of  the 
Langevin  Equation  (2.12),  (2.13)  by  considering  propagation  of  \hz  field  statistics  directly,  using 
the  concept  above  of  the  equivalent  deterministic  medium  (EDM),  as  embodied  in  the  operator 
Qi(d\  etc.,  cf.  (2.34).  For  example,  since  =  Qi(d).  we  may  write  (2.4),  taking  the  average 

of  both  sides,  as 


<(£(0)  -  0)  a  >  =  <aH>  =  och  ,  (2.40a) 

or 

(£(°)  -  0i(d))  <ot>  =  aH ,  (2.40b) 

which  is  called  a  Transport  Equation,  or  here  a  propagation  equation  of  the  mean  (total)  field,  <a>. 
From  (2.37b)  et  seq.,  it  is  evident  that  (2.40b)  is  an  integro-differential  equation  in  <a>.  Second- 
and  higher-order  generalizations,  for  the  propagation  of  <ai0t2>,  etc.,  follow  in  similar  fashion. 
(For  a  treatment  of  the  propagation  of  field  statistics,  see  [13],  [14],  Chapter  14,  for  instance.) 

2.7  An  Exact  Solution  for  a  and  <a> 

Under  certain  conditions  it  is  possible  to  obtain  exact  solutions  for  the  inhomogeneous  (and 
total)  field  ai,  a.  Physically,  this  can  be  done  for  infinite  media  where  boundary  interactions  can 
be  ignored,  which  may  be  geometrically  possible,  so  that  only  volume  effects  occur.  We  consider 
the  case  of  deterministic  inhomogeneities.  It  is  then  required  that  the  (deterministic)  inhomogeneity 
operator  Q  takes  the  form 

•e 

&R,t  IR’,t%  =  JdR'  /  dt,A(R-R’,t-t')(  )R*,r 
(->  — 

=  f  Je'ik,R+st  (  )mAo(M)  (2.41) 

J  (2ti)3  BJ  2m 
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where  Ao(k,s)  =  F~j(F  t{  A(r,t)},  and  the  kernel,  A,  of  @  (230),  is  thus  of  convolutional  type  in 
space  and  time,  but  #a8(R-R')8(t-t'),  and  is  accordingly  nonlocal.  The  result  is 


a(R,t)  =  FkF~l  (Oo)  =  J 

Ck 


dk  rGoT(t,s)yo(k,s)^at,R-|-sl 
(2ti)3  B  J  1  ~  N0(k,s) 


Nq  =  ~  Do,«®  Aq 


ds 

2tci 


(2.42) 


where  %i00  is  given  by  (2.25a,b).  Here  Ck  is  an  appropriate  contour  in  the  complex  k-plane  to 
ensure  that  the  radiation  condition  ("spatial  causality")  is  obeyed  [IS].  Bri  denotes  a  two-sided 
Laplace  transform,  or  Bromwich  contour  (-«>+id,  «»+id),  d  £  0,  which  enables  us  to  handle 
transient  signals  and  signals  of  finite  duration  (d  >  0),  as  well  as  the  steady-state  cases  (d  =  0), 
which  latter  reduce  to  the  temporal  Fourier  transform  situation,  cf.  Sec.  2.2.  The  evaluation  of 
(2.42)  is  now  reduced  to  quadratures,  which  may  usually  be  carried  out  analytically,  at  least  in  an 
appropriate  series  form. 

For  unbounded  random  media  the  same  approach  may  be  followed.  We  reduce  the 
governing  Langevian  equation  to  a  Dyson  equation,  as  above,  cf.  (2.41),  where  now  — >  ^2 

and  a  -»  <a>.  Equation  (2.42)  now  applies  if  we  replace  a  by  <a>,  N0(k,s)  by  (k,s)  = 
-%(k,s)oaq<jd)(k,s),  where  q&d](k,s)  s  F  R1Ft{Qi(d*)»  with  the  kernel  of  Qi(d)  having  the  form 
q^OR-R'.t-t'),  cf.  (2.41)  above.  For  mixed  deterministic  and  random  media,  cf.  (2.1),  a  linear 
combination  of  operators  Qi(d)  with  the  above  properties  may  be  expected. 

2.8  The  Received  Field,  X(t) 

So  far,  we  have  been  considering  the  (total)  field  a(R,t),  at  some  point  P(R,t),  cf.  Fig. 

2.1.  For  subsequent  signal  processing  this  field  must  be  sensed  by  an  appropriate  receiving 
aperture,  or  array,  designated  by  the  (linear)  filter  operator  ft,  viz. 

X  =  fta,  or  X  =  {Xi, ....  Xj}  =  fta,  (j  =  m,n);  J  =  MN,  (2.43) 

where  (2. 12)  now  represents  the  received  wave,  in  continuous  or  discretely  sampled  form  in  space 
and  time.  Beam-forming  is  achieved  by  summing  over  the  spatial  inputs,  e.g.,  X  =  {EmXjsmj,). 
See  Section  IYB  of  [4]. 

In  fact,  whether  we  are  concerned  with  remote-sensing  of  the  medium  or  detection  and 
classification,  we  must  inevitably  include  the  receiving  sensor  (s),  embodied  in  ft.  Accordingly,  in 
such  cases  we  must  replace  a,  <a>,  etc.,  by  X,  <X>,  etc.,  in  the  above  and  subsequent  analyses. 
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This  is  formally  straightforward,  but  will  introduce  additional  computational  requirements.  We 
remark  that  from  the  viewpoint  of  the  transmitter-receiver  systems  of  Fig.  1. 1  the  fields  we  are 
dealing  with  here  are  all  signal-dependent ,  i.e.,  they  are  generated  by  a  designated  active  source 
(Sin,  Fig.  1.1)  and  can  produce  random,  signal-dependent  noise  (reverberation  or  incoherent 
scatter)  in  the  receiver.  This  is  in  contrast  to  ambient  noise  fields,  whose  mechanisms  are  quite 
independent  of  the  designated  source  Sin.  For  the  purpose  of  the  present  treatment  we  regard  the 
always-present  ambient  noise  as  essentially  negligible;  but  see  Sec.  IV,  A,  [4]:  the  effects  iof 
ambient  noise  sources  may  be  formally  including  if  we  now  replace  the  source  function  ~Gj  by 
-Gt-<ja,  when  Ga  is  (also)  a  distributed  source  mechanism. 

3.  Diagram  Methods  and  Approximations 

An  important  aid  in  integrating  and  approximating  the  formal  solutions  of  the  field 
equations  here,  and  in  achieving  numerical  results  eventually,  are  diagram  methods  ([3],  [12],  Part 
I,  [16]).  These  enable  us,  along  with  the  relevant  physics,  to  select  the  dominant  contributions  and 
omit  the  unimportant  ones,  and  often  to  choose  approximations  which  considerably  reduce  the 
ensuing  calculations.  We  illustrate  the  approach  with  first-order  results  in  the  case  of  random 
media. 

Note  that  for  deterministic  media  these  r  ^uns  apply  as  well  if  we  drop  the  various  averages 
over  the  Langevin  ensemble ,  whirh  new  contains  only  one  member  field  representation,  (Xq. 

3.1  Basic  Diagrams 

We  begin  with  the  following  ensemble  diagrams,  for  the  ensemble  FOR,  FOS,  PTSS,  etc., 
and  FD  (Feynman  diagram)  equivalents  for  the  field  ensemble  [a]  at  any  allowed  point  P(R,t): 


FOR: 

1  >'  ■■  O  ( >  >  I 

v- «, — '  6  > — 1 

cf.  Fig.  2.1 

(3.1»'' 

FD: 

1=3 

*  W/////Jk  ♦ 

H 

__  X  [ - 1  ; 

!  //' 

Eqs.  (2.10,  2.12) 

(3.1b) 

(») 

*  >H  * 

».iw  ; 

FOS: 

I 

;  1 - 1  •-*©  ■H.l-tyl 

M  .  -6.  ] 

Eq.  (2.12) 

(3.1c) 

PTSS: 

•  CZ) 

■[i»  *— o  ■ 

N 

(3.1d) 

Here  we  have  — >  (s  Kfc*)  representing  th e  feedforward  operator,  while  <— 0  (=  {Q})  denotes  the 
ensemble  feedback  or  random  scattering  operator,  Q.  The  ensemble  Feynman  diagram  equivalent 
(FD)  of  the  FOR,  (3.1a),  is  precisely  (3.1b),  while  (3.1c)  gives  the  corresponding  ensemble  FOS, 
and  (3.  Id)  the  PTSS,  which  is  readily  obtained  by  iterating  the  FD,  (3.1b),  or  as  the  unaveraged 

23 


TD  8589 


form  of  (2.33).  The  (deterministic)  first-order  Dyson  equation,  Eq.  (2.35)  is  now  easily  written: 
First-order  Dyson  Equation'. 

H 

<(X>  =  Oh 

where  now  we  denote  the  average  by  <  >,  applied  to  the  corresponding  ensemble  symbol.  Our 
general  diagram  conventions  are  to  represent  ensemble  quantities  as  "open"  forms,  e.g.,  I  I , 
4-0,  o,  etc.,  and  average  and  otherwise  deterministic  quantities  by  solid  forms,  e.g.,  4— •,  •,  •, 
with  field  averages  denoted  by  brackets,  <  >,  e.g.,  <==>  =  <a>,  etc.  The  aim,  of  course,  is 
consistency  of  notation  and  simplicity  with  a  minimal  evolution  of  complexity. 

Diagrams  may  be  manipulated  in  the  same  fashion  as  are  operator  equations  [cf.  (2.35), 
etc.]:  one  observes  the  same  positioning  rules  and  conditions  of  inversion  [cf.  (2.13)],  e.g.,  for 
multiplication  and  "division."  Addition  and  subtraction  are  also  directly  equivalent  to  the  addition 
and  subtraction  of  operators  and  operator  derived  (algebraic)  quantities,  as  indicated  in  the 
examples  below.  Finally,  transforms  of  diagrams  are  the  corresponding  diagrams  of  the 
transforms.  This  permits  an  explicit  diagrammatic  representation  of  the  various  wave  number- 
frequency  (WNF)  spectra  of  the  scattered  field,  particularly  the  amplitude  and  intensity  spectra  in 
operator  form. 

Accordingly,  the  corresponding  diagrams  for  the  equivalent  medium  [cf.  Sec.  2.5]  become 
[cf.  (3.1)]: 


Equiv.  Med.  FOR: 
(2.35) 


Equiv.  Med.  FOS: 
(2.35) 


Equiv.  Med.  FD: 

First-Order  Dyson  Equations : 
(2.35) 


(3.3a) 


(3.3b) 


(3.3c) 
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PTSS: 

(2.13) 


FOS  +  PTSS: 
(2.12),  (2.13) 


***’  •  O-lM;  •»<$>;  (3.3d) 

h.«j>  <(n.(})'2>>  n  t, 

•  Cl  ♦  T/fZ/V/A 


<*=>  •  i-1 — •  ci  .'mm.  n 

1-  H  4  t  H  V  / 

4 


Here^”'1  ^\n^2)  denotes  the  n-fold  average;  •  s  <Q>;  ■  <^(2)>  .  <Q0£)O>, 

k  <^0^oo^n>,  etc.  Equations  (3.3e)  follow  from  the  iteradon  of  the  FD  (3.3c).  In  the 
same  way  we  find  that  the  EDiM,  cf.  (2.34),  (2.35)  as  developed  in  (2.37)  with  (2.37a,b),  can  be 
expressed  diagrammatically  as 


;u) 


(3.4a) 

(3.4b) 


this  last  in  terms  of  r\x  (s  «=0=>  ),  where,  for  example,  <n0>  =  <*=0=  *Co>  etc. 


32  Diagram  Approximations 

A  major  problem,  as  always  for  strong-scatter  situations,  is  the  approximate  evaluation  of 
the  FOS  or  its  PTSS  equivalent,  since  exact  analytic  solutions  are  possible  only  in  special  instances 
[cf.  Sec.  2.7].  Another  important  problem  is  our  ability  to  truncate  the  PTSS  in  comparatively 
weak-signal  regimes:  for  example,  what  is  the  effect  of  stopping  at  a  given  term  in  the  PTSS? 
What  is  the  overall  contribution  of  the  remaining  terms?  What  criteria  of  truncation  are  reasonable, 
etc.?  In  any  event,  various  approximations  are  in  order. 

We  consider  briefly  here  two  useful  classes  of  approximation:  I.  Series  modification ,  and 
II.  Truncation.  In  the  former,  the  character  of  the  equivalent  deterministic  medium,  or  "mass" 
operator  is  modified.  In  the  latter,  the  PTSS  is  truncated:  a  finite  number  of  interaction  orders  is 
retained,  while  the  higher-orders  are  discarded.  For  example,  let  us  consider  the  former  case  (I) 
first: 


A.  Series  Modification:  First  Order  Dyson  Equation 

As  applied  to  the  various  forms  of  the  Dyson  Equation  [cf.  (2.35),  (3.2)],  the  form  of  the 
equivalent  deterministic  scattering  operators  (EDSO),  Q,(d),  §[d2\ .  .  .,  <$$m,  cf.  (2.34),  is 
modified  in  some  fashion,  suggested  by  the  physics  of  the  problem.  This  EDSO  is  altered,  usually 
by  truncation.,  while  the  infinite  operator  series  implied  by  (1  -T)^-1),  etc.,  cf.  (2.35),  etc., 

25 


TD  8589 


remains.  Let  us  examine,  for  example,  the  first-order  Dyson  equation  (2.35)  for  the  mean  field 
«x>  and  postulate  a  first-order  independent  structure  for  & 

/.  First-Order  Independence: 


:  =  •  «... 


’n » 


(3.5) 


and  from  (3.5)  in  (2.37a,b)  and  then  in  (2.35)  we  see  at  once  that  the  EDjSO,  Q|d),  is  directly 


(3.6) 


where  0id)  reduces  to  a  single  term,  so  that  (2.35)  becomes  at  once 

«x>  =  an  + <^>  <a>  =  (f  -  Kt*,  <Q>)-]  an 

with  the  following  equivalent  (deterministic)  diagrams  [cf.  (3.2),  (3.3)]: 

♦“■>  •  ; — ^ Mi  •  "j  ~ • 

1  "  H  1  '  H  H 

Thus,  in  the  case  of  purely  volume  scatter  (<Q>  =  0).  Q(d)  =  0,  and  <oo  =  an,  here,  which 
clearly  shows  that  the  postulate  (3.5)  is  much  too  restrictive  to  provide  meaningful  results  in  most 
cases. 

In  a  similar  way,  we  may  employ  a  more  structured,  second-order  assumption  on  the 
statistical  character  of  0-  This  is  the  second-order  "independent"  or  (generalized)  Bourret  or  Bi¬ 
local  approximation,  defined  by 


(3.7) 


(3.7a) 


26 


TD  8589 


II.  Second-Order  Independence 


^  (3.8a) 

: -  ^  ^  ^(n)  •  (3.8b) 


For  the  odd-order  case  we  can  also  set,  alternatively, 

n  (n-l)/2 

<m>  ~  <^>  n<C>i(2>  ;  n  =  odd  :  =  •  ^  ^  ...  ^(n) .  (3.8c) 

j  i 

Consequently,  the  EDiSO,  Q(d)),  (2.37),  (2.37a, b),  becomes  here 

+  <—  At(1)  -»  «-  A^  — >  +  ...: 

£h(d)  |B  =  <0>  +  [  <$&>.&>  -  <0>^u<0>]  -  [<Q>tooo<Q&*&>  -  <0>^oo<0>^l»<C>]  + ... 

(3.9a) 

or 

Q[d)  Ib  =  <Q>  +  ( -  <Q>Aoo<0>]  -  [<Oi0u^i0u<€>  -  <§>&*,<&>&„<&>]  + ... 

(3.9b) 


<m> 

j 

n/2 

=  n<0P>  ;  n  =  even 

j=i 

(n-l)/2 

=  n<Qi(2»<^>  ; n  =  °dd 

i 

(Oi=0  =  t)  • 

In  diagram  form  these  are: 


dll 


C 


(3.9c) 


Now,  unlike  the  purely  independent  case  (3.5),  (3.6), ...  the  series  for  Q|d)  does  not  terminate,  so 
that  in  this  second-order  situation  we  are  generally  forced  to  truncate  the  series  to  obtain 
manageable  approximations. 

Note  from  (2.37b)  that  when  <$  >  =  0,  which  represents  purely  volume  scatter  effects 
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without  boundaries  or  interfaces,  the  bi-local  or  Bourret  approximation  (3.%)  for  reduces  to  a 

single  term  -. 


<&>  sisQ; _ _ 

Ibo^b  =  <£M>  <-  A(»);  A  < »  -  0.  m  a  2:  O  . 


(3.10) 


for  this  first-order  equivalent  deterministic  medium.  [As  we  have  just  seen  (cf.  (3.9a,b)),  with 
interfaces  the  Bourret  approximation  (3.8)  yields  a  non-terminating  series.]  This  compact  result 
suggests  the  utility  of  the  bi-local  approximation  in  applications,  wherever  it  appears  reasonable  to 
describe  the  principal  scattering  effects  of  the  medium  through  second-order  spatial  (and  temporal) 
correlations  at  most. 

The  various  mean  fields  <ot>  are  diagrammed  here  in  (3.3)  when  Q[d))  iB  is  given  either  by 
(3.9)  or  by  (3. 10)  in  the  unbounded  (unlayered)  cases.  For  the  latter  we  have  explicitly 


0: 


I  Bourret 


i 


mm  t 
"  tn[d)i,  • 


(3.11) 


for  the  associated  FOS  and  PTSS.  Unlike  the  strictly  independent  case  I  above,  here  <a>  *  an 
and  <oo  now  represents  a  meaningful  expression  for  the  mean  scattered  field  when  there  are  no 
interfaces,  cf.  remarks  following  Eq.  (3.7a).  Equation  (3.1 1)  is  the  generalization  of  Tatarskii’s 
Eq.  (26),  ^61,  [17],  based  on  the  "classical"  Helmholtz  equation  (2.6),  with  time-independent  £  = 
e(R,t),  to  arbitrary  (linear)  media  and  input  signals.  The  corresponding  analytic  solution,  similarly 
generalizing  Eq.  (26a)  of  Tatarskii  (ibid.)  is  given  here  by  (2.42),  modified  for  random  media. 
Sec.  2.7. 


B.  Series  Truncation 

In  the  cases  above  the  "mass  operator"  for  EDSO,  etc.,  is  approximated  either  as  an 
infinite  or  terminated  series  [cf.  (3.9),  (3.6),  (3.10)],  thus  altering  the  form  but  not  the  fact  that  the 
resulting  approximate  solution  is  still  an  infinite,  though  sub-series  of  the  original  PTSS.  On  the 
other  hand,  truncation  replaces  the  exact  (or  previously  approximated)  infinite  operator  series  with 
a  finite  operator  series.  Of  these  latter  there  are  two  principal  types:  I:  Taylor's  series 
approximations,  which  yield  analytic  forms  directly,  and  H:  Born  approximations,  which  constitute 
a  hierarchy  of  truncated  operator  series. 

Let  us  consider  I  first: 
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/.  Taylor  Series  Approximations: 

Although  the  exact  solution  <a>  of,  say,  the  first-order  Dyson  equation  (2.35),  (3.2), 
requires  an  infinite  PTSS,  cf.  (2.25),  we  can  reduce  the  evaluation  of  «x>  to  a  purely  algebraic 
solution  by  means  of  a  Taylor's  series  expansion  with  truncation.  Thus,  if  we  expand  <a>: 

<a(R',t')>  =  <a(R,t)>  +  (R'-R)  •  Vr  <a(R,t)>  +  (t'-t)  ^  <a(R,t)>  +  ...  ,  (3.12) 

and  keep  only  the  leading  term,  we  have  directly  the  first-order  approximate  form 

[&oA(d)  <«>]  =  <a(R,t)>  lCLA(d)l  (3.13a) 

where  AoJ^i(d)l  is  simply  an  (analytic)  function  of  (R,t),  viz.: 

r\[d)l  =  ]0u4(d)l  =  /  Moo(R,t  IR'.t')  dR'dt'  J  <$  {d)(R\t*  IR",t")l  dR"dt"  s  n  }d)(R,t).  (3.13b) 


Applying  (3.13a,b)  to  (2.35)  gives  the  approximate  analytic  result 


<aO)(R,t)>  = 


«H(R.t) 

1  -  n{d)(R,t) 


(3.14) 


where  the  exact  or  various  approximate  expressions  for  the  equivalent  deterministic  medium 
operator  Q(d)  may  be  employed,  cf.  ((2.37a,b),  (3.6),  (3.10)).  Thus,  for  the  fully  independent  (I) 
and  bi-local  approximations  (II),  <£>>  =  0,  we  have  specifically  (and  exactly)  here: 

nUd  =  &-<$>!;  <0>  =  0.  (3.15) 

Qearly,  there  may  be  some  points  (foci)  or  lines  (caustics)  in  space  where  <a^)>  —» oo,  when  the 
denominator  of  (3.14)  vanishes.  The  reality  of  these  singularities  depends  on  how  well  (3.14) 
approximates  <a>  itself:  the  true  foci,  etc.  may  not  be  infinite,  but  may  be  located  close  to  the 
singularities  of  the  approximation.  Thus,  (3.14)  may  serve  as  a  guide  as  to  potentially  real  foci, 
etc. 

Finally,  the  error  in  using  (3.14),  say,  is  generally  difficult  to  assess,  because  the  true 
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of  «x>  is  not  usually  available.  However,  in  the  cases  where  the  conditions  for  (2.41)  are  obeyed 
we  can  use  the  exact  result  (2.42)  to  determine  the  error  vis-i-vis  using  (3.14)  for  <a>. 


II.  Born  Approximations: 

Unlike  the  above  procedures  A,  8,  which  essentially  approximate  the  "mass  operator," 

Q(d),  with  an  infinite  or  truncated  series,  cf.  (3.9)  or  (3.10),  we  may  truncate  the  PTSS  [(2.13), 
(3.3d),  etc.]  directly.  The  resulting  finite  series  is  called  a  Bom  approximation  series,  whose  order 
depends  on  the  stage  at  which  truncation  is  applied.  Thus,  for  die  first-order  Dyson  equation 
(3.3c)  for  the  mean  field  <a>  we  can  construct  the  following  hierarchy  of  approximations  and 
associated  diagrams:  mu),: 

n  =  0:  ++  4  — ^ ►©  ‘  V  (3.16a) 


n-1: 

first  Born  approximation 

n  =  f 

second  Bom  approximation  I 

etc. 

The  n1*1  order  FOR  is  simply 
(FOR)n: 

nth-order  Born  approximation 


4^  i  [!♦—*•]  VM/A-.  (ros)l:<»>i 

i  [!♦  — *•  ♦  E23.  trOJ)2: 

2  M 

«x>  =  Z<(^-^)(n)>aH ; 

ne?c.  etc.  (n£3). 


n»Q 


«■! 


■  <a>:  ■ 

I  <(l'1..0)(B)>»M. 
n»0 


(3.16b) 

(3.16c) 


(3.16d) 


Unlike  the  approaches  above  which  use  modificaticms  of  the  "mass  operator,"  Q^,  etc.,  but  where 
the  resulting  PTSS  remain  infinite,  or  where  equivalently,  the  corresponding  FOR  and  FOS  [cf. 
(3.1 1)]  contain  closed  feedback  loops,  all  Bom  approximations  contain  only  feed-forward  loops, 
as  indicated  in  (3.16d).  This  is  the  immediate  consequence,  of  course,  of  the  truncation  of  the 
direct  PTSS  (3.3d). 


Examples :  Ocean  Volume  and  Surface  Channels:  The  Received  Field  X 

As  examples,  we  noted  that  both  the  important  cases  of  weak  volume  scatter  in  the  ocean 
and  ocean  surface  scatter  are  well  approximated  by  at  most  the  scattering  element  interacting  with 
the  homogeneous  field  cch,  the  former  because  multiple  scatter  effects  are  usually  quite  ignorable, 
the  latter  from  the  physical  geometry  of  the  wave  surface,  which  likewise  discourages  multiple 
scatter  except  within  a  typical  wave  crest-to-crest  domain  (except  possibly  at  small  grazing  angles), 
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cf.  Fig.  2.1,  (n  =  1).  Similar  considerations  usually  apply  for  ocean  bottoms,  as  well,  cf.  Fig. 

2.1.  We  shall,  however,  consider  only  the  modes  M  =  S,  V  here,  observing  also  that  0(30-40  db) 
they  are  effectively  independent  (see  [4],  Ref.  21,  thesis]):  there  is  negligible  coupling  in  the  ocean 
between  surface  (bottom)  and  volume  scatter.) 

Accordingly,  a  first-Bom  approximation  suffices,  and  we  write  for  the  received  field,  X, 
cf.  (2.43), 

Qm 

X  =  &a  =  &  ( — >0  +  -m-0-»Q)  =  ft(ctH  +  K^MOtH)  -*  :  &  «-0:  £>M 
A  S  or  V  :  A 

M  =  S,  V.  ©  =  -Gt  (3.17a) 

The  mean  wave  (X)  is,  with  <  >  =  •> 

<X>  =  ft<a>  =  ft(  — »©  +  — >t— »  0)  =  {  ft  (an  +  ftAc^NP^H)  }  •  (3.17b) 

XH  +  Xi 

The  second  moment  of  X,  (3.17a),  is  found  to  be,  on  averaging  over  the  phases  (<J>)  of  the  input 
signal. 


«XiX2»(j>  =  <ft  ift2<aia2»(j> 


=  +  I^1^2  *  |<0iO)0iO)>s  +  <qi(1)qlu>  +  X<qi(k)qik)> |aHiaH2] ^ ^ 

*•  k=2  J  (3.18a) 


*  =  0 

k  -  1 

(k  £  2) 

{&M‘\ 

+  — -  oCb 

* 

+  ►W)GZZ3), 

11  22 

11  22 

1 1  22  H,j 

_  "homog."  . 

*  unseat.  ' 

"classical":  coherent 

+  incoherent 

multiple-scatter  incoher. 

FOM:  coh. 

+  incoherent 

(0) 

where 


O  h  q  =  0  -  <0>R  =  *-0-0  =  L  jv(0)  =  X(<v(k)>R  +  Av(k)> 


k=l 
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&°)  =  <&>R  =  O,  <(^(°)>S  =  <Q>  =  •  =  <0>R 
•  \  s<  >S:space av. x 0 

=  <  qi  q2>  etc.  =  <  >  =  <  >r,s  (3.19) 

1  2 

Note  that  the  (total)  average,  <  >,  as  in  [30],  Appendix  A.l  et.  seq.  involves  both  an  average  over 
radiation  events  and  an  average  over  random  spatial  (S)  and  parameter  values  (6),  embodied  in 
<  >s,  cf.  (3.14)  and,  particularly,  [30],  Appendix  A.2. 

It  is  the  decomposition  principle,  embodied  in  Qm.  (M  =  S,  V,  B),  ([30],  (A.2-5)  et  seq.), 
which  allows  us  to  resolve  the  inhomogeneity  operator  into  sets  of  distinct  and  (statistically) 
independent  entities  as  specifically  exhibited  in  (3.18a)  and  (3.18b).  Moreover,  and  most 
important,  the  entities  (k  =  0,  k  =  1,  k  2: 2)  in  (3.18a)  and  (3.18b)  above  and  in  (3.17),  with  the 
help  of  (3.19)  have  an  explicit  physical  interpretation.  Thus,  the  terms  k  =  0  contain  both  the 
coherent  and  incoherent  radiation  contributions,  from  all  orders  (k  £  1)  of  radiation  interactions, 
e.g.,  single  and  multiple  scattering  (k  1 2),  for  example.  Consequently,  if  vOO  represents  the 
density  of  (illuminated)  k-coupled  scattering  elements  (k  >  1),  then  <vOO>r  is  the  coherent 
radiative  contribution,  where  <  >r  is  the  average  over  the  radiation  events  associated  with  the 
ensemble  of  potential  (re-)radiating  sources.  In  a  similar  way 


AvOO  =  vOO  -  <vOO>r  (3.20a) 

is  the  fluctuation  in  the  density  of  k-coupled  radiating  elements,  and  is  always  associated  with 
incoherent  radiation.  Accordingly,  we  have 

•O  M 

v(0) s  ^(<v(1c)>r  +  AvOO)  =  svi <v(°)>r  +  Av(°)  (3.20b) 

k=l  k=l 

showing  the  basically  coherent  and  incoherent  elements  of  the  scattering  or  source  region 
containing  these  various  k  £  1  types  of  elements.  The  (average)  density  of  coherent  radiators  is 
<v(°)>r  =  L<vOO>r,  since  <Av('£)>r  =  0.  The  number  of  "radiation  events"  occurring  in  a  small 
region  dA  of  the  illuminated  or  emitting  domain  A  is  dN,  so  that  v  =  dN/dA  =  v(°).  A  taxonomy 
and  interpretation  for  the  dN,  similar  to  (3.20a),  (3.20b),  can  be  constructed:  one  now  has 

M  *0 

dN(°)  s  £(<dN(k)>R  +  dnOO)  =  ^dNto.  (3.20c) 

k=l  k=l 

The  relations  (3.21)  et  seq.  obviously  apply  if  one  replaces  v  etc.  by  dN. 
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Since  the  medium  in  question  is  linear,  the  scattering  or  inhomogeneity  operator  Q  is  a 
linear  functional  of  the  vOO,  i.e.,  of  v(°),  e.g.,  $  =  -£(v(°)),  thus  containing  all  types  of  radiation 
interaction  effects.  Moments  of  v(°),  or  more  generally,  of  ^m,  which  appear  in  (3.17)  and 
(3.18),  above,  are  readily  indicated:  we  have 

<Q>  =  £{<v(°)>) 

where 


<v(°)>  -  «v(°)>r>s.  (3.21) 

This  is  different  from  zero,  in  scattering,  if  specular  reflection  geometry  (both  on  a  surface  an  din 
a  volume)  is  available.  The  second-order  moments  are  instructive: 

<Qi&>  =  £{<v{0)v|0)>)  (3.21a) 

where  now 

•o-pv^0)>  =  «Av1(0)Av^0)>r>s  +  «v  ](0)>r<v^0)>R>S-  (3.21b) 

with 

Av(0)  h  ^AvOO. 
k=l 

Equation  (3.21b)  contains  both  incoherent  (~  Av(°))  and  coherent  terms  (~  <Av(°)>r),  of  all  orders 
(k  2: 1)  of  coupled  scattering  or  radiating  elements.  The  associated  covariance  of  Q  is 

Kq=  £{<vfM0)>-  <v{0)><v£0)>}  (3.21c) 


which  embodies  wholly  incoherent  radiation. 

As  in  all  these  approximate  situations,  the  evaluation  of  error  remains  a  major,  very  difficult 
problem,  because  the  exact  results  are  generally  unknown.  We  must  rely,  instead,  on  our  physical 
intuition,  aided  by  the  (essentially)  local  physics  governing  the  inhomogeneity  operator  0  in 
question.  This  is  one  reason  why  the  physical  "anatomization"  of  the  so-far  canonical  0  is 
critically  important  in  specific  cases. 

33  An  Estimation  Procedure  for  the  Mean  Operator 

When  the  elements  (-Gt,  $L)  of  the  "input-output"  structure  of  Fig.  2.2  are  available, 
estimates  of  the  "mass  operator"  may  be  made  to  any  degree  of  accuracy  from  the  empirical  or 
a  priori  mean  field  <a>,  in  a  controlled  way,  in  principle,  with  determinable  error.  This  is  a  form 
of  "system  identification"  problem,  typical  of  control  theory  [24].  Here,  in  essence,  one 
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Fig.  3.1  Feedback-operational  method  for  obtaining  the  (first-order)  equivalent  medium  (integral) 

operator 


introduces  a  "black-box"  or  system  with  a  controllable  input,  observes  the  resulting  output(s),  and 
attempts  to  infer  the  operations  of  the  system  which  relates  the  two.  A  key  feature  here  is  that  the 
system  in  question  is  known  to  be  linear,  so  that  unique  relations  are  established. 

For  (2.34)  to  be  a  physically  (vis-it-vis  mathematically)  useful  device,  we  must  be  able  to 
indicate  a  physical  procedure  for  determining  Qid)-  This  may  be  done  according  to  the  scheme 
shown  here  in  Fig.  3.1.  This  scheme  produces  a  converging  series  of  estimates,  Qj^t’  f°r  0{d). 
by  iteration,  when  [<a>,  -Gj )  are  provided  experimentally,  and  when  iCt*,  is  either  empirically  or 
analytically  given.  The  basic  approach  is  conceptually  quite  simple:  with  a  starting  estimate  Oisi-i 
of  0{d),  one  runs  through  the  second  loop,  using  the  (experimentally)  observed  «x>,  the  given 
source,  -Gj,  and  the  known  global  descriptions  of  the  homogeneous  component  of  the  medium 
(Aoo  here,  cf.  Sec.  2.2),  to  get  an  <a>est-  This  is  then  subtracted  from  the  observed  <a>  from  the 
first-loop,  which  in  turn  leads  to  an  adjustment  of  a  second  comparison,  with  a  further 
reduction  in  the  error  A<a>  =  <a>  -  <a>est  and  A($[d)  *  Qi(d)  -  0{^t-  Various  optimum  and 
suboptimum  schemes  from  control  theory  (e.g.,  Wiener  filtering,  gradient-climbing,  etc.,  now, 
however,  extended  to  four  dimensions)  may  be  used  to  drive  the  errors  A<ot>,  A^(d)  essentially  to 
zero,  so  that  ->  <3  jd).  This  technique  may  also  suggest  acceptable  approximations  in 
experimental  situations.  In  fact,  approximating  here  is  analogous  to  choosing  some 
approximation  to  Q[d)>  (2.37),  like  Ao,  or  Aj,  cf.  (2.37  a,b,c)  et  seq. 

In  all  cases  the  particular  physics  of  the  problem  will  dictate  the  procedural  details,  the 
approximations  (usually  required),  and  the  class  of  "hardware"  to  be  employed,  e.g., 
minicomputer,  minisupercomputer,  or  supercomputer.  Some  form  of  the  latter  appears  necessary 

34 


TD  8589 


in  most  cases,  cf.  Sec.  4  ff.,  since  these  problems  are  generally  on  the  scale  of  long-range  weather 
prediction.  In  fact,  these  problems  may  be  considered  as  a  form  of  "underwater  weather" 
prediction,  especially  when  the  medium  (and  interface)  behavior  is  properly  included.  It  is  hoped 
that  the  preceding  material  may  provide  at  least  a  stimulus  to  the  much-needed  and  as  yet 
nonexistent  computational  solutions  to  this  important  class  of  problems. 

4.  Concluding  Remarks:  Analytic  Steps  to  the  "Interface  Problem" 

As  we  have  noted  at  the  beginning  of  this  paper,  a  principal  purpose  of  this  study  is  to 
bring  the  formal  analytic-operational  treatment  to  the  point  where  the  computational  problems  of 
implementing  these  formal  results  numerically  can  be  addressed.  We  call  this  the  "Interface 
Problem,"  where  the  interface  lies  between  the  general  mathematical  formulation,  outlined  here, 
and  software  design  and  hardware  performance. 

The  various  steps  which  can  bring  the  analytical  treatment  to  this  interface  (the  dotted  line) 
are  summarized  in  Table  4.1. 

Our  philosophy  here  has  been  to  start  with  as  general  a  "top-down"  approach  as  possible, 
including  those  physical  features  which  can  ultimately  influence  practical  solutions.  Among  such 
features  are  beam  patterns,  signal  waveform,  and,  of  course,  geometry,  as  well  as  the  relevant 
structures  of  the  medium  inhomogeneities.  Proceeding  with  particular  applications,  we  usually 
find  that  approximations  and  simplifications  can  be  made  which  can  significantly  reduce  the  very 
heavy  computational  burdens  implied  by  the  operational  forms  of  the  desired  "solutions."  The 
latter,  of  course,  are  both  deterministic  and  stochastic,  and  a  full  treatment  will  need  to  include  both 
these  components.  The  present  problems  are  on  the  scale  of  weather  prediction  problems.  In  fact, 
they  are  themselves  models  of  "underwater  weather."  With  the  recent  advances  in  computer 
capabilities  (and  reduced  costs)  we  may  expect  these  problems  to  be  within  range  of  numerical 
solution:  The  present  effort  has  been  presented,  in  part,  to  stimulate  their  numerical  investigation, 
as  well  as  the  achievement  of  further  analytic  solutions. 
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Table  4.1 


Step  1  Deterministic  Inhomogeneities 


1 .  Obtain  the  PDE  and  b.c.'s  +  i.c.'s ;  Eq. 
(2.1)  et  seq.;  =  "Dynamical  Equation";  all 
representations  are  equivalent  here 


Stg_ 


1. 

la. 


2.  Obtain  the  infinite-medium  Green's  function 
goo.  cf.  (2.24),  (2.26).  (2.27) 

3 .  Operational  Solutions :  These  are  complete: 
PTSS,  FOS  from  the  FOR,  cf.  (2.13),  Fig. 
2.2„  Sec.  2.3,  2.4  for  the 
inhomogeneity  operator,  cf.  (2.5b),  (2.6a). 
(2.9) 


2. 

3. 

3a. 


3b. 


4.  Approximate  0.  local  and  distributed. 

Use  diagram  methods ,  Sec.  3.2 

[0  is  usually  local,  cf.  (2.5b),  (2.6a),  (2.9)] 

5.  Analytic  solutions  ( for  volumes  only),  cf. 
Sec.  2.7,  Eq.  (2.42) 


6.  Other  Approximations:  Ray  tracing,  normal  6. 
modes  (time-independent)  [21];  Q(R  I,..) — 
time-harmonic  solutions;  Q(t  I...}—  space 
harmonic  sol.;  & R  ,t  I...) — sol.  not  harmonic. 
Medium  non-reciprocal,  so  Generalized 
Huygens  Principal  not  available.  Path 
Integrals,  Super-Eikonal  Methods  [22],  [23] 

7.  Computational  Solutions:  special  cases  only  7. 
(no  boundaries)  [7],  [10].  General  software: 

to  be  estimated,  approx,  and  exact. _ 


I  Random  Inhomogeneities _ 

Obtain  the  PDE,  b.c.'s  +  i.c.'s;  Eq. 

(2.1)  et  seq.,  for  a  typical  representation 
Construct  the  Lange vin  Equation,  Eq. 
(2.4)  et  seq. 

Same 

Solutions  are  the  statistics  of  the  field  ot, 
e.g.,  «x>,  <aict2>,  etc.,  and  pdfs 
wi(a),  W2(ai,a2), ... 

Construct  Dyson's  equation  (first- 
order),  Sec.  2.5  for  Equivalent 
Deterministic  Media  (EDM),  sol.  <a> 
Higher  order:  £)^,  for  <ai<X2>,  etc. 
PrSS;  FOS  forms  (2.35);  Some  analytic 
solutions.  Sec.  2.7 
Approximate  etc..  Sec.  2.5; 
Diagram  methods  assist;  Sec.  3,  Sec. 

3 2,  is  non-local,  all  higher-order 

frd)’s  are  also  non-local 
Analytic  solutions  (volumes  only),  cf. 
Sec.  2.7,  Eq.  (2.42)  modified,  2nd  ft 
1st  order.  Higher-order  by  similar  way, 
but  very  complicated. 

Same,  for  (EDM),  (EDM)2,  via 
etc.  (See  [12]  for  details.)  Realization 
of  diagram  approximations 


Computational  Solutions :  generally 
unknown.  To  be  constructed,  approx. 
and  exact  =  Major  Problem. _ 
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Glossary  of  Principal  Symbols 


а,  a(R,t) 

Ot,  Or 
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Vc 

V2 

e,  e(R,t) 

M(d) 

FD 

FOR 

FOS 

FRO 

lFi 

77-1 

б,  G 

goo 

»  *  K 

*x>  *y>  •z 
k,  k 

t(0) 

L 

lOuC^oo 
N0,  Ng? 

PTSS 

<3(d).  Os.V.B.... 

R 

ft 

Sin 

X,  X(t) 

%,oo.  To  ,oo,  ^  0,°° 

t 


=  total  field,  at  P(R,t) 

=  beam  patterns 
=  homogeneous  field 

=  Bromwich  contour,  (two-sided  Laplace  transform) 

=  mean  s  ound  speed 
=  sound  speed  gradient 
=  Laplarian  (in  rectangular  coordinates) 

=  (function  of)  index  of  refraction 
=  field  renormalization  operators  (FRO),  (EDFRO) 

=  Feynman  diagram 
=  feedback  operational  representation 
=  feedback  operational  solution 
=  Field  Renormalization  Operator 
=  confluent  hypergeometric  function 
=  Fourier  transforms 
=  imput  source  operators  and  functions 
=  Green's  function  of  an  unbounded  medium 
=  unit  vectors 
=  wave  numbers 

=  differential  operator  for  homogeneous  part  of  medium 
=  linear  operator,  cf.  Eq.  (2.22) 

=  integral  Green's  function  operator 
=  field  renormalization  functions 
=  perturbation  theoretical  series  solution 
=  inhomogeneity  operators 
=  position  vector  in  rectangular  coordinates 
=  receiving  array  operator 
=  input  signal 

=  (complex)  angular  frequency 
=  medium  operator 

=  receiver  array  outputs  to  signal  processor 
=  Fourier  and  Laplace  transform  of  g« 

=  unit  operator 
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Prof.  H.  R.  Raemer  1 

IDA-CRD 

Dr.  D.  I.  Lieberman  1 

Dr.  J.  DeLucia  1 

IDA-SRC 

Dr.  P.  Schneck  1 

Dr.  H.  Freitag  1 

IDA,  Alexandria 

Dr.  R.  D.  Turner  1 

LAB  FOR  PHYSICAL  SCIENCES 

Dr.  T.  Beahn  1 

UNIV  OF  ARIZONA 

Prof.  J.  R.  Wait  1 

UNIV  OF  CAMBRIDGE.  Cambridge.  England 

Dr.  B.  J.  Uscinski  1 

Dr.  B.  Kerman  1 

Prof.  Longuet-Higgins  1 

MULTIFLOW  COMPUTER.  INC 

Dr.  C.  Pangali  1 

DEFENCE  RESEARCH  ESTABLISHMENT 

Dr.  M.  F.  Chapman  1 

DR.  D.  MIDDLETON.  New  York  25 

DR.  A.  O.  SYKES.  Vienna.  VA  1 

DR.  M.  SCHULKIN.  Potomac.  MD  1 

DR.  D.  HOLLIDAY,  Sacramento.  CA  1 

DR.  D.  P.  PETERSEN  1 

DR.  W.  A.  VON  WINKLE  1 
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